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Abstract 

Using an idea going back to Madelung we construct global 
in time solutions to the transport equation corresponding to the 
asymptotic solution of the Kolmogorov- Feller equation describing 
a system with diffusion, potential and jump terms. To do that 
we use the construction of a generalized delta -shock solution of 
the continuity equation for a discontinuous velocity field. We also 
discuss corresponding problem of asymptotic solution construction 
(Maslov tunnel asymptotics). 
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Introduction 

The goal of the present paper is to present a new approach to the con- 
struction of singular (i.e., containing the Dirac 5- function as a summand) 
solutions to the continuity equation and to show how these solutions can 
be used to construct the global in time solution of the Cauchy problem 
for Kolmogorov-Feller- type equations with diffusion, potential and jump 
terms. The relation between the solutions of the continuity equation and 
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the system consisting of the Hamilton-Jacobi equation plus the trans- 
port equation has been well studied before in the case of a smooth action 
functional. 

The velocity field u is determined as the family of velocities of points 
on the projections of the trajectories of the Hamiltonian system corre- 
sponding to the Hamilton-Jacobi equation. As it was mentioned by E. 
Madclung [13], in this velocity field, the squared solution of the transport 
equation (denoted by /rho) satisfies the continuity equation 

Pt + iV,up) + ap = (0.1) 

with some additional term ap, which is defined below (a is equal to if 
the Hamiltonian is formally self-adjoint). The main obstacle to the ex- 
tension of this correspondence globally in time is the fact that in general 
the solution of the Hamilton-Jacobi equation are smooth only locally in 
time. The loss of smoothness is equivalent to the appearance of singu- 
larities of the velocity field mentioned above. Till recent times there 
was no method for constructing formulas for solutions of the continuity 
equation for a discontinuous velocity field. It is clear that the continu- 
ity equation has the divergent form and this very important property 
allows precisely to introduce the concept of a global solution in spite of 
singularities in the velocity field. The divergent form as itself does not 
play any role in Madelung's approach but it is very important for our 
construction because we deal with singular solutions. 

In the present paper we generalize Madelung's approach to the case 
in which the singular support of the velocity field is a stratified manifold 
transversal to the velocity field trajectories. This holds, for example, in 
the case where the space is one-dimensional under the condition that, 
for any t G [0,T], the velocity field singular support is a discrete set 
without limit points. 

It is of interest to note that the generalization of Madelung's idea, 
first suggested for the Schrodinger equation, also makes sense for equa- 
tions from a different class. 

Namely, let A" be the Lagrangian manifold corresponding to the 
solution and let tt: A" ^ be the projection mapping. 

All the points of 7r~^(x) G A" make contributions to the construction 
of an asymptotic oscillating solution (of the type of the WKB-solution 
(p{x, t) exp(fS'(x, t)/h)). 

But if we consider nonoscillating solutions of the form ^p{x, t) exp(— s(X, r)/if),| 
just as for equations of heat conduction type with a small parameter h 
at the second-order derivative, see Sec. 2 and [8], then only the point 
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from the set 7 G tt ^{x) at which 

S {gamma, t) = min S{'y,t) 

gammaE.7r~ ^{x) 

makes a contribution to the solution, hence the points of the Lagrangian 
manifold making contributions to the construction of a nonoscillating 
solution form a surface with (shock wave-type) jumps, which results in 
discontinuities of VS and hence in discontinuities of the velocity field 
for the transport equation. 

The last equation does not have a divergent form in contrast to 
the continuity equation, which arises in the Madclung construction. So 
Madelung's idea was initially adopted to the global constructions involv- 
ing generalized solutions to the continuity equation. 

In Section 2, we describe the definition and construction of these 
types of solutions. The conclusion is that these solutions can be con- 
structed by means of characteristics. This allows us partly to change the 
direction of time and to solve some inverse problems. We will discuss 
this problem in detail in subsequent papers. 

1 Generalized solutions 

of the continuity equation 

Here we follow the approach developed in [1] , where the solution of the 
continuity equation is understood in the sense of an integral identity, 
which, in turn, follows from the fact that relation (0.1) can be understood 
in the sense of the distributional space P'(M"t^). The first step in this 
way has been done in [14], see also [15], [16], where the approach based 
on smooth approximations of the solutions was used. 

We specially note that the integral identities in [1] can be derived 
without using the construction of nonconservative products [2, 4] of the 
nonsmooth and generalized functions (or measure solutions [5]), and the 
value of the velocity on the discontinuity lines (surfaces) is not given 
a priori but is calculated. In the case considered in [1], the integral 
identities exactly coincide in form with the identities derived using the 
construction of a nonconservative product (measure solutions) in the 
situation described at the end of above introduction, which we shall 
now make more precise. 

First, we consider an n — 1-dimensional surface 74 moving in M^, 
which is determined by the equation 

It = {x;t = ip{x)}, 
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where G C''^(]R"'), and VV' 7^ in the domain in M" where we work. 

This is equivalent to determining a surface by an equation of the 
form 

S{x,t) = 

{S G in both variables, S{x, t) = 0, Vx,tS\s=o 7^ 0) under the condi- 
tion that 

We remain that the situation with ^ = 0, can also be covered 
by making the change of variables x\ = Xi — Cit with appropriately 
chosen Cj, i = l...n, solving the problem with the moving surface and 
then returning to the original variables. Possible generalizations are 
considered later in this section. 

Next, we assume that ( belongs to C^(M" x R^). Then, by defini- 
tion, 

{S{t - ^{x)), C(x, t)) = I ^{x)) dx, 

where 6 is the Dirac delta function and (, ) is the distributional pairing 
(with respect to the variable t G R\ and x G M"). 

Let S{t — il){x)) be applied to the test function 77 G Co°(M"), then 



{6{t - ip{x)),r]{x)) = r)du^, 
J It 



where du) is the Leray form [6] on the surface {t = ^^{x)} such that 
—dtpduj^ = dxi . . . dxn- 

One can show that (see [1], [6]) 

{6{t - ^{x))Mx)) = 1^ ^^da. 

First, we assume that the solution p to Eq. (0.1) has the form 

p{x, t) = R{x, t) + e{x)6{t - tp{x)), (1.1) 

where R{x, t) is a piecewise smooth function with possible discontinuity 
at {t = ■4){x)}: 

R = Ro{x,t) + H{t - i/j{x))Ri{x,t), 

e G C(M"') and has a compact support, ip & C"^ and Vip 7^ for a; G 
suppe, and H{z) is the Heaviside function. 
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It is clear that the term 

e{x)5'{t-tp{x)) 

appears in (0.1) if we differentiate the distribution S{t — tp{x)) with 
respect to t. Hence it is necessary to have in (0.1) 

(V,pu) = —e{x)S'{t — + smoother summands, 

since VS{t — ip) = —Vip6'{t — tp). Then we must have 

— 1 — TTo ~ V') + smoother summands. 

I V'i/'r 

Now we formulate an integral identity, defining a generalized solution 
to the continuity equation. 

We set r = {{x,t)\t = V'(^)}; this is an n-dimensional surface in 



u{x,t) = uo{x,t) + H{t — il^)ui{x,t), 



r X Mj^. Let 



where tp is the same function as before, and uo,ui G 

Let us consider Eq. (0.1) in the sense of distributions. For all C{x,t) G 
C^(M" X R\_), C(x, 0) = 0, we have 

+ (V, pu),c) = -{P, Ct) - {pu, VC). 

Substituting the singular terms for p and pu calculated above, we come 
to the following definition. 

Definition 1.1 A function p{x,t) determined by relation (1.1) is 
called a generalized 5-shock wave type solution to (0.1) on the surface 
{t = V'(^)} if the integral identity holds 



{RCt + {uR, VC) + aRQ dx dt 
for all test functions C{x,t) G C^{W x R\), C{x,0) = 0, ^ = 

(f^>v) + |v^I|. 

We have also the relation 
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We note that the vector n± is orthogonal to the vector (VV',— 1), 
which is the normal to the surface F, i.e., lies in the plane tangent 
to r. 

We can give a geometric definition of the field . The trajectories 
of this vector field are curves lying on the surface F, and they are or- 
thogonal to all sections of this surface produced by the planes t = const. 
Furthermore, it is clear that the expression is an absolute value of 
the normal velocity of a point on 7^, i.e., on the cross-section of F by 
the plane t = const, and the expression • =^ Vn is the vector of 
normal velocity of a point on 7^. Thus, we nave another representation: 



where Vn = 7r*(u„), u„ is the normal velocity of a point on 7^, and vr* is 
induced by the projection mapping tt : F — ^ E^. 

It follows from the latter definition that the following relations must 
hold: 

Rt + (V, Ru) + aRC = 0, for all points {x, t) ^ F, 



{[R]-\Vi^\[Run]) + \^—j ^=0' for all points (x,t)GF, 

(1.3) 

The last relation can be rewritten in the form 

]CE + -^E= [RunWl - [R], (1.4) 
an 

where E = e/|VV'|, the factor /C = (V, y^^) = div v {v is the normal on 
the surface {t = tp{x)}) and, as is known, is the mean curvature of the 
cross-section of the surface F by the plane t = const, ^ = ( |^^, V). 
Now we assume that there are two surfaces 

Ti = {{x,t)\t = il)i{x)} 

in X M]j_, z = 1, 2, whose intersection is a smooth surface 

7 = {(x,t);(t = Vi)n(t = V2)} 

belonging to the third surface F(3) = {{x,t);t = iIj^{x)}. Further, we 
assume that the surface F(3) is a continuation of the surfaces F^*) in 

the following sense. We let ny denote curves on the surfaces Fj and 
we assume that each point (x,t) on the surface 7 is associated with 
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a graph consisting of the trajectories n^f^ and n^^ entering {x,i) and 

the trajectory n!f^ leaving this point (i.e., the trajectories n^J^ fiber the 
surface r*-*^). We also assume that the surface (stratified manifold) Fu = 
T^i-j U r(-2) U r(3) consists of points belonging to these graphs. Next, we 
assume that u{x,t) is a piecewise smooth vector field whose trajectories 
enter Ty. 

Definition 1.2. Let 

3 

u{x, t) = uo{x, + ^ H{t - 'ipi)uii{x, t), 
1=1 

where ip is the same function as before, and uo,uii G C{W^ x M5|_). The 
function p{x, t) determined by the relation 

3 

p{x,t) = R{x,t) + '^ei{x)6{t - ijji{x)), 

i=l 

where R{x,t) G C^{W x R\) \ {Ijrf^}, is called a generalized (5-shock 
wave type solution to (1.2) corresponding to the stratified manifold Fy 
if the integral identity 

roc p 

I 1 {RCt + {uR, VC) + aRC) dx dt 
Jo Jr" 

holds for all test functions C{x,t) G C^(R" x R\), Cix,0) = 0, 

(Ffe>V) + |V^.||. 

As above this relation implies the first equation from (1) outside Fy, 
equations of the type of the second equation in (1) on strata F(j) and 
the Kirchhoff type relation on 7: 

(ei + 62)1- =63!-. (1.6) 

Now we consider the case with codimF > 1. First, we note that the 
second integral in (1.2) can be written as 



dn 



L 



' -ax, ^) = // ( v) + ^) c(x, t) dx. 



dn± 



We note that if the surface F is determined by the equation S{x,t) = 
rather than by the simpler equation {t = il^{x)} presented at the 
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beginning of this section, then 



\vs\ \vs\ \vs\^ 

and, of course, the new vector field = {Vn, V) + ^ remains tangent 
to r. 

Therefore, in this more general case, using this new vector Vn, we 
can again rewrite the integral identity from Definition 1.1 as 

[ [ {RCt + {uR,VC)+aRC)dxdt+ [ e({Vn,V) + ^)C{x,t)dx = 0. 
Jo jRri Jr ^ Cr/ 

(1.7) 

This form of integral identity can easily be generalized to the case 
in which F is a smooth surface in M"+^ of codimension > 1. 

In this case, instead of Vn, we can use a vector v that is transversal 
to r and such that the field {v, V) + ^ is tangent to T. We note that the 
vector V is uniquely determined by this condition, which can be treated 
as "the calculation of the velocity value on the discontinuity" from the 
viewpoint of [5] and [7]. 

Moreover, in this case, the expression for p does not contain the 
Heaviside function, and it is assumed that the trajectories of the field u 
are smooth, nonsingular outside F, and transversal to F at each point 
of F. In this case, the function p has the form 

p{x,t) = R{x,t) + e{x)S{r), 

where R G Ci(M'*+^ \F), e G C^{T), and the function d{r) is determined 

by 

(<5(r),c(x,t)) = / (io, 



where uj is the Leray form on F. If F = {Si{x,t) = 0(1- • •r\Sk{x,t) = 0}, 
k G [l,n], then lo is determined by the relation, see [6], p. 274, 

dtdxi dxn = dSi dSkLO. 

In this case, we assume that the functions are sufficiently smooth 
(for example, C^(R"' x M^j.)) and their differentials on F are linearly 
independent. 

Moreover, we can assume that the inequality 

j_ V{Si,...,Sn) 

T){t,Xi, . . .,Xn-l) 
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holds. This inequahty is an analog of St 7^ at the beginning of this 
section and allows us to write u in the form 



OJ = J ^dxk dxn- 

The integral identity, an analog of (1.7), has the form 

J {RCt + {uR,VO + aRC) dxdt + J e(^iv,V) + ■^^C{x,t)u = 0. 

Integrating the latter relation by parts, we obtain equations for deter- 
mining the functions e and R similarly to (1.4). 

Now we assume that the singular support of the velocity field is the 
stratified manifold |J Fj with smooth strata Ft of codimensions > 1. 

We also assume that the velocity field trajectories are transversal to 
(Jr and are entering trajectories. 

Then the general solution of Eq. (0.1) has the form 

p{x,t) = R{x,t) + 

where R{x, t) is a function smooth outside |J Fj, are functions defined 
on the strata Fj, and the sum is taken over all strata Fj. 

The integral identities determining such a generalized solution have 
the form 



poo r 

I / {RCt + {uR, VC) + aRC) dxdt 

Jr" 



+1: 



ei 



V^,V) + ■|-^]ax,t) 



iOi = 0. (1.8) 



This implies that, outside IJ^i, the function R satisfies the continuity 
equation 

Rt + (V, uR) +aR = 0, 

and, on the strata Fj for nj = 1, equations of the form (1.4) hold, 
which contain the values of R brought to Fj along the trajectories. For 
ni = n — k, k > 1, on the strata F;, we have the equations 

d 

Q^eifii + (V, vieiin) = Fini, (1.9) 

where m is the density of the measure ui with respect to the measure 
on F; which is left-invariant with respect to the field ^ + (i7;,V), and 
Fi is defined by the following construction. Denote a e-neighborhood of 
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Ti by and denote its boundary by dVf. Let us consider the integral 
appearing after integration by parts: 

JdTf 

where is the normal component of the velocity u on dTf, ujf is the 
Leray measure on dTf, C is a test function. Passing to the limit as e ^ 
we obtain 



It is well known that outside |J Fj the function R{x,t) can be calcu- 
lated using the famous Cauchy formula 

R(x.t) = po(x,t)\^\~^ eM- f adt') (1.10) 
Dxq Jo 

where po ia a constant along the trajectories of the field u outside |J Ti, 
is the jacobian of the mapping corresponding to the shift along 
the trajectories of u and the integral under exponent is calculating along 
the trajectories of the field u. 

This formula implies that the limit as e ^ of the above integral 
exists. 

We note that it follows from the above that the function R is deter- 
mined independently of the values of Vi on the strata under the condition 
that the field trajectories enter (jFj. 

In conclusion, we consider the case where the coefficient a has a 
singular support on [jTi, i.e., 

a = f{u). 

In this case, we set 

ap = ap + ^ f{vi)ei5{Ti). 

where d = f{u) outside IJ Fj. We note that such a choice of the definition 
of the term ap is not unique in this case. But, first, it is consistent with 
the common concept of measure solutions (sec [3], [5]) and, second, it is 
of no importance for the construction of the solution outside \jTi for 
the case in which the trajectories u enter (jTj. 
In this case, identity (1.8) takes the form 

/ / {RCt + iuR,VO + f{u)RC)dxdt 

Jo JR" 

+ 1^/, ei[((^;i,V) + - + /(t;.))c(a;,t)]cu, = 0, (1.11) 
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and Eq. (1.7) can be rewritten in the form 



^(ezw) + {V^vieini) + f{vi) = Fim. 



(1.12) 



All the afore said gives the following statement. 
Theorem 1. Let that the following conditions he satisfied for t G [0, T], 



(1) yjVi is a stratifies manifold with smooth strata Fj/ 

(2) the trajectories of the field u are smooth outside (jTi, enter (jFj 
and do not intersect outside (jFj/ 

(3) equations (1-12) are solvable on the strata Fj/ 

(4) the Kirchhoff laws are satisfied on the intersections of strata Fj . 
Then there exist a general solution to the continuity equation (0.1) 

with a = f{u) in the sense of the integral identity (1.11). 

2 The Maslov tunnel asymptotics 

We recall that the asymptotic solutions of a general Cauchy problem for 
an equation with pure imaginary characteristics was first constructed by 
V. P. Maslov [8]. In the present paper, we consider only the following 
Cauchy problem 



-h^ + pil, -h— ]u = Q, u{x, t, h)\t=o = e-^''W/V°(^), (2.1) 



where P{x,$,) is the (smooth) symbol of the Kolmogorov-Feller operator 
[9], 5*0 > is a smooth function, (p^ £ , h +0 is a small parameter 
characterizing the frequency and the amplitude of jumps of the Markov 
stochastic process having having transition probability given by P(x, ^). 
To be more precise, we can have in the mind the following form of 



where A{x) is positive definite smooth matrix, V{x) is a smooth function 
and its second derivatives are assumed to be uniformly bounded, and 
Pl{x, dv) is a family of bounded measures smooth with respect to x in the 
sense that the functions x ^ iJ,{x, B) are smooth for all measurable sets 
B. The symbol P{x,^) can also depend on t, we will be more precise 
later on. 



T > 0: 



1 




P{x,0: 
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Locally in t, an asymptotic solution of problem (2.1) can be con- 
structed according to the scheme of the WKB method, see [8]: the 
solution is constructed in the form 



u = e 



oo 

-S{x,t)/h 



in the sense of asymptotic series. In this case, for the functions S{x, t) 
and ipo{x,t) we obtain the following problems: 

^ + ^'(x,^J =0, S{x,t)\t=o = Soix), (2.2) 



d^o r.( dS\ ^ \ d'^P d'^S , , 



^3 



the Jacobian 



As is known, the solution of problem (2.2) is constructed using the 
solutions of the Hamiltonian system assumed to exist and to be smooth 

x = V^P{x,p), x\t=o = xo, (2.4) 

p = -VxP{x,p), p|t=o = VS'o(xo). 

This solution is smooth on the support of (po{x,t) for all t such that 

Dx/Dxq 7^ for xo G supp ip^{x). We let denote the 

translation mapping along the trajectories of the Hamiltonian system 
(2.4). 

We recall that the plot 

= {x = xo,p = VSoixo)} 

is the initial Lagrangian manifold corresponding to Eq. (2.2), and A" = 
5(^Aq is the Lagrangian manifold corresponding to Eq. (2.2) at time t. 
Let TT : A" — > be the projection of AJ? on M^, which is assumed to 
be proper. For this property to hold, it is sufficient to assume that the 
trajectories of the system (2.4) do not go to infinity during a finite time. 
The point a G AJ* is said to be essential if 

S{a,t) = min S{P,t) 

/367r-l(Q;) 
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and nonessential otherwise. Here S is the action on A" determined by 
the formula 



Jo 



pdx — H dt, 



where the integral is calculated along the trajectories of the system (2.4) 
the projection of its origin being xo = P- As is known 

S{x,t) = S{Tr-^x,t) 

at regular points where the projection vr is bijective. 

The global in time asymptotic solution of problem (2.1) is given by 
the Maslov tunnel canonical operator. 

To define this operator, following [8, 10] we introduce the set of 
essential points \Jyit C A". This set is closed because the projection vr 
is proper, i.e. that for all x the set of p such that {x,p) G A", Tr{x,p) = x 
is finite. 

Suppose that the open domains Uj C A" form a locally finite covering 
of the set {jju- If the set Uj consists of regular points, then we set 



= (>Sj(x,t)/h 



ipoj{x,t) 



(2.5) 



where 



ipoj{x,t) = ll^Oj{x,t) 



Dxq 



Dx 



1/2 



ipOj{x,t) being the solution of the equation 

^ + {P^ix, VSj), V^oj) - \ ^(^, V5,)^o,- = 0, 



(2.6) 



which exists and is smooth whenever 



Dx/Dxq 



^ 0. The solution 



Uj in the domain containing essential (nonregular) points (at which dn 
is degenerate) is given in the following way: the canonical change of 
variables is performed so that the nonregular points become regular, 
then we determine a fragment of the solution in new coordinates by 
formula (2.5) and return to the old variables, applying the "quantum" 
inverse canonical transformation to the solution obtained in the new 
coordinates. 

The Hamiltonian determining this canonical transformation has the 
form 



1 



o-kPk, 
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where cri, . . . cr„ = const > 0. 

The canonical transformation to the new variables is given by the 
translation by the time —1 along the trajectories of the Hamiltonian 
H„. One can prove (see [8], [10]) that the family of sets a for which the 
change of variables takes a nonregular point into a regular is not empty. 

Next, the solution near the essential point is determined by the re- 
lation 

Uj = eT^^'^Uj, (2.7) 
where Uj is given by formula (2.5) in the new variables and 




On the intersections of singular charts (containing singular points) 
and nonsingular charts (without singular points), we must match Sj and 
ipQj. This can be done by applying the Laplace method to the integral 
whose kernel is a fundamental solution for the operator —h^+H^. This 
integral appears if we write the right-hand side of (2.7) in detail. In this 
case, since the solution is real, the Maslov index which is well-known 
[8] to appear in hyperbolic problems does not appear. The complete 
representation of the solution of problem (2.1) is obtained by summing 
functions of the type (2.5) and (2.7) over all the domains Uj, for more 
detail, see [8], [10]. Here we only say the corresponding sum is (locally) 
finite because we assumed that the projection vr is proper. 

The asymptotics thus constructed is justified, i.e., the proximity be- 
tween the exact and asymptotic solutions of the Cauchy problem (2.1) 
is proved [8, 9]. More precisely it is proved that at the points of the set 
7r(U 7it) where the projection tt is bijective the following estimate holds: 

u{x, t, h) — Uj = 0{h) 

In the preceding Section we noted that the values of the solution 
of the continuity equation at nonregular points are independent of the 
values of the solution on the singularity support (of course, the inverse 
influence takes place) by the condition that the velocity field trajectories 
enters the singular support. 

In the case of the canonical operator whose construction has briefiy 
been described above, the relation between the solutions at essential 
and nonessential point is also unilateral, namely, the essential points 
are "bypassed" using (2.7), but the values of the functions Uj on the 
singularity support do not determine the values at the regular points. 
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It is clear that the values of the functions Sj and tpoj at regular 
points arc defined by characteristics via the initial data. The trajectories 
of the Hamiltonian system also come to points on A" such that the 
projection mapping at these points is singular. At these points, we 
cannot define V'oj by characteristics directly and we must use an auxiliary 
construction (see (2.7)) which allows us to determine the values of an 
asymptotic solution at the projections of singular points onto x-space. 
In this auxiliary construction, the values of Sj and ipoj in a neighborhood 
of the singular points on A" are used to determine the values of Uj at 
the projections of singular points. Thus, there is a similarity between 
this and the preceding sections: we define the values of an asymptotic 
solution by characteristics outside its singular support and then define 
the values of the asymptotic solution at the singular points using the 
already defined values at regular points in our auxiliary construction. 

Now we note that the function S{x,t) is such that 

S{x,t)\u, = Sji7r-\a),t) 

is globally determined and continuous at points of the domain 7r(|J ju) C 
M^. We denote this set by (jFj and assume that this is a stratified 
manifold with smooth strata Tu of different codimensions. We note that, 
for example, if the inequality V{Si{x,t) — Sj{x,t)) ^ holds while we 
pass from one branch A"n(J to another, then the set 7r{{Si — Sj) = 0} 
generates a smooth stratum of codimension 1. In the one-dimensional 
case, all strata are points or curves on the (x, t)-plane (under the above 
assumptions about the singularities being discrete). 

Now we consider the equation for tp^j . We denote this function by p 
and then obtain 

^ + {V,up)+ap = 0, (2.8) 

where u{x, t) = V^P(x, V5) and a = — tr §^{x, V5). 
If the condition 

Hess^P(x,^) > 

is satisfied, then it follows from the implicit function theorem that 
VS{x,t) = F{x,u{x,t)), where F{x,u) is a smooth function and 

a = f{x,u), 

where f{x, z) is again a smooth function. 

Thus, we have proved the following theorem. 

Theorem 2. Suppose that the following conditions are satisfied for t G 

[o,r], r>0; 
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(1) There exists a smooth solution of the Hamiltonian system (2.4). 

(2) The singularities of the velocity field u = V^P(x, VS') form a 
stratified manifold with smooth strata and Hess^ P{x, > 0- 

(3) There exists a generalized solution p of the Cauchy problem for 
Eq. (2.8) in the sense of the integral identity (1.9). 

Then at the points ofTT{\Jjit) where the projection tt is bijective, the 
asymptotic solution of the Cauchy problem (2.1) has the form 

u = exp(-5(x, t)/h){^ + 0{h)). 

This theorem is a global in time analog of the corresponding Made- 
lung observation [13] about local solutions of Schroedinger type equa- 
tions. 



3 Particular cases 

The theorem stated in the previous section requires that some assump- 
tions are satisfied. The most restrictive is the item 3 in theorem 2, that 
is the existence of the global generalized solution to the continuity equa- 
tion. Under the above made assumptions it is possible to construct this 
solution using characteristics, but only in the case where the structure 
of the singular support of u is not changing in time-all sections of the 
stratified manifold introduced above by planes t = const smoothly de- 
pend on t. A more complicate situation arises when the singularities of 
the velocity field change their structure. In this case the problem of the 
construction of a global in time generalized solution to the continuity 
equation has not been solved yet. The obstacle is that in this case usu- 
ally one has no global in time expression for the velocity field u. In turn 
this does not allow to apply formula (1.10) to construct global solution 
to the continuity equation. In the multy-dimensional case as far as we 
know there is only one result concerning to shock wave generation [14] 
which allows to construct global in time approximations of the shock 
wave formation process. But this is slight different from the construc- 
tion that we needs here. In the one dimensional case the situation is 
better and we have all needed formulas. 

We begin with the spatially homogeneous case. Here the problem is 
equivalent to the one of constructing a formula for a global solution to 
conservation law equation 

Here P{—h-^,t) is the same operator as in (2.1) but assumed to 
be independent of x with the symbol P{^,t) and v = dS/dx. The 
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velocity field u in this case is P^{v, t). In [11] a construction of the global 
solution to the continuity equation where the velocity field is given by 
the solution of the equation (3.1) was given. Because the set of singular 
points is discrete by our assumptions, without loss of generality one can 
consider the case were only one point of singularity appears. Denote the 
corresponding (smooth) initial condition by uq, the instant where the 
singularity appears by t* and the point of singularity by x*. 

The first step of construction suggested in [11], [12] is that we 
change uq in a small neighborhood of of origin Xq of the trajectory 
coming to x* when t = t*. We denote this new part of initial data 
ui{xo) for xq G {xq — I3,Xq + (3), P —>■ and assume 

ep-'^ ^0, e ^ 0. (3.2) 

We define the function ui = Ui{xo, t) as a solution of implicit equa- 
tion 

P^{ui,t) = -K{t)xo + b{t), (3.3) 

The latter equation is solvable under the condition Hess^P(x,^) > 0, 
formulated above. 

The functions K(t) and b{t) are defined from the condition of conti- 
nuity of the characteristics flow, i.e 

ui{xq - I3,t) = uo{xq - 13, t), ui{xq + /?,*) = uo{xQ + P,t) 

It is easy to check that this choice of ui provides that the Jacobian 
|Z)a;/£)xo| is identically equal to for t = t* and xq G [xq — /?, Xq + 
/?). Here we remove from usual topological concept of general position 
considering the situation of identical equality that can be destroyed by 
small perturbation. But this construction follows from the algebraic 
concept and allows to present the solution of (3.1) in the form of linear 
combination of Heaviside functions (see [11]). 

The second step of our construction of an approximation is a mod- 
ification of the definition of characteristics. We set 

x = {l-B)P^{ui{xo,t),t) + Bc, xoe{x*o-l3,x*o + P) (3.4) 

and 

X = P^{uo,t), 
where xq does not belong to {xq — P,Xq + f3), 

P{v{x{xl + t), t)) - P{v{x{xl - /?, t), t)) 
v{x{xl +P,t),t)- v{x{x^ - 13, t),t) 
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The initial data for (3.4) are as follows: 



= xo + Ae, £ > 0. 

=0 

The function B in (3.4) has the form B = B{{t - t*)/e) and B(z) is 
smooth, monotone and increasing from to 1 for z G (—00,00). Simi- 
larly to [11], [12] one can prove that there exist an ^ = const such that 
the Jacobian \Dx/Dxo\ calculated using the above introduced charac- 
teristics is not equal to zero, but it is of order 0{e) when t > t* + O{0) 
when xq G {xq — (3,Xq + j3). Using the velocity field generated by x we 
can construct global in time (smooth) solution of the continuity equa- 
tion in the form (1.10). After that, passing to the limit as e ^ we will 
obtain the generalized solution of the continuity equation in the sense 
of definition from Sec.l just like it was done in [12]. 

Spatially inhomogeneous one dimensional case. 

We will follow the scheme introduced above. The case under con- 
sideration can be treated in the same way as the previous one with 
modifications. Firstly, we will assume that the symbol P = P(x,4) does 
not depend on t. In this case this assumption (which means that the 
mapping g^p is invertible w.r.t. time) will be used to construct the in- 
sertion to initial data. In the previous case we did it using the implicit 
function theorem, see (3.3). 

Let Aq be a smooth nonsingular (w.r.t the projection tt) curve in 
the {x,p) space, which is a Lagrangian manifold corresponding to initial 
data for our problem. We consider the Lagrangian manifold A|* = g^p Aj 
and assume that there is only one point singular with respect to the pro- 
jection onto X-axis and its projection is x* . Let (3 be the same as above. 
Let us set t\ = t* + (3. Because of the assumption that P'^^ is positive, 
we have that for t = t\ the Lagrangian manifold Aj, has two parts which 
contain essential points and these parts form a shock wave type curve 
with the jump at the point x\ where Sif,ft{x\,t\) = Sright{x\,t\). We 
connect these parts by a vertical line and thus obtain a new Lagrangian 
manifold, which is a piccewise smooth continuous curve with two angle 
points (ends of the vertical part, the distance between them of order 
/?). We denote this manifold by Aj, and apply the mapping gj!'^ for 
sufficiently small ti to this manifold. This mapping obviously exists and 
is a diffcomorpfism because our Hamiltonian P does not depend on t. 
We consider the obtained manifold g^^^ Aj* as the new Lagrangian man- 
ifold corresponding to our problem for t = t\— ti changing the manifold 
Aj* by gip*^ A^, . As it was said above the latter manifold is piecewise 
smooth curve with two angle points and all points of the curve outside of 
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the part between these angle points are regular. Moreover there exist a 
sufficiently small ti such that the part of the curve between these angle 
points contains only regular points-these statements are a consequence 
of the positivity of P", its stationarity and the possibility to choose ti 
small enough (and independent on e). 

Denote the projections of the mentioned above angle points on the 
manifold gp^^A^* to the x-axis by ai < a2- We note that |ai — 02] is of 
order /3. 

Like in the previous example we introduce the new characteristics 
system 

x = {l- B)Pl{x{xo,t),p{xo, t)) + Be, (3.5) 
p = -{l- B)P^{x{xo,t),p{xo,t), xo e (01,02), 

and 

x = p^(x,p), (3.6) 
P = -P'xix,p) 
when Xo does not belong to (01,02). We have set 

^ ^ P{v{x{a2.,f^;P{a-2,f))) - P{x{ai,t),p{ai,t)) 
p{x{a2, t),t) - p{x{ai,t),t) 

The initial data for (3.5), (3.6) are as follows: 

x\^^Q = xo + Ae, 
p\t=o =Poixo), 

where (xq, po{xo))=gp^'^ Aj:, . The expression on the right hand side of 
(3.7) is the direct analog of the well known Rankine-Hugoniot expression 
for the velocity of the shock propagation. In the case under consideration 
it is the velocity of the point x on x-axis, where Sieft{x, t) = Sright{x, t). 

By the assumption we have only one singular point if we are consider- 
ing the family of manifolds A|, t G [0, t*\. We also have by construction 
that the Jacobian J = Dx/Dxq calculated using the solutions of the 
system (3.5) is not equal to zero. More precisely we have 

lim J = H{t* - i) Jo, 

where Jq is the Jacobian calculated using the solutions of (3.5) for S = 
( Jq = when t = t* by construction) and 

J > H{t* -t)Jo + Ce, 
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where C = const > 0. This statement directly fohows from (3.5) if we 
take the properties of the function B into account. This means that the 
velocity field, generated by the projections of the solution of the system 
(3.5), (3.6) onto the x-axis has nonintersecting trajectories for £ > 0. 
Thus we can use it to construct solutions of the continuity equation. It 
remains to note that just like in [12] it is easy to check that the limits of 
these solutions will satisfy to the integral identities introduced in Sec. 1 
as the definition of generalized solutions to continuity equation. 
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